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We summarize the renormalization-group properties of transverse-momentum dependent (TMD) 
parton distribution functions (PDF)s arguing that in the light-cone gauge the overlapping ultra- 
violet and rapidity divergences cannot be solely controlled by (dimensional) regularization, but 
necessitate their renormalization. In doing so, we show that at the one-loop order this additional 
divergence entails an anomalous dimension which can be attributed to a cusp in the gauge con- 
tour at light-cone infinity. Then, we present a recent analysis of TMD PDFs which incorporates 
in the gauge links the Pauli term ~ F^ v [Yn,Yv]- This generalized treatment of gauge invariance 
is shown to be justified, in the sense that it does not modify the behavior of the leading-twist 
contribution, though it contributes to the anomalous dimension of that of twist-three. An im- 
portant consequence of the inclusion of the spin-dependent Pauli term is the appearance of a 
constant phase — the same for the leading twist-two and subleading distribution functions — that 
ensues from the interaction of the Pauli term in the transverse gauge link with the gauge field 
accompanying the fermion. Remarkably, this phase has opposite sign for the Drell-Yan process 
as compared to the semi-inclusive DIS. 
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1. Introduction and Theoretical Framework 

One of the problems inherent in the definition of hadronic observables is how to ensure gauge 
invariance. This problem arises because correlators are nonlocal quantities which contain local 
operators that transform differently under gauge transformations, hence entailing a dependence 
on the gauge adopted. Clearly, physical quantities should not depend on the choice of the gauge 
we choose to work — this should merely be a matter of (calculational) convenience. To render 
integrated parton distribution functions (PDF)s gauge invariant, it is sufficient to insert into their 
definition a Wilson line — a gauge link — between the two Heisenberg quark operators that renders 
their product gauge invariant |JlJ]. From the point of view of renormalizability, this operation intro- 
duces additional contributions to the anomalous dimension of the PDFs. These contributions stem 
from the local obstructions of the gauge contours: endpoints, cusps, and self-crossing points (see 
[^] for a technical exposition and references). It is to be emphasized that although the gauge link is 
nonlocal, no explicit path dependence is introduced, e.g., on the gauge-contour length. Actually, to 
ensure the gauge invariance of the PDFs it is even sufficient to use a straight lightlike line, because 
integrated PDFs are defined on the light cone and the only contribution from the gauge link to the 
anomalous dimension of the PDF comes from its endpoints (see, for instance, [^] and references 
cited therein). Hence, one has for the integrated PDF of a quark i in a quark a 

yi/ fl W = ^/^-^«-<p|^(r,o ± )y f [r,o-|«1w(o-,o ± )|p> ) (i.i) 



where 



[£-,0-\tf] = ^exp 



ig ( dzW &z- ,0x)t a 

J0-[V] 



(1.2) 



is a path-ordered gauge link (Wilson line) in the lightlike direction from to £ along the contour H>. 

One may insert a complete set of states and split the gauge link [£~,0 _ ] into two gauge links 
connecting the points and £ through °°. This is mathematically sound, provided the junction 
(hidden at infinity) of the two involved contours is smooth, i.e., entails only a trivial renormalization 
of the junction point so that the validity of the algebraic identity [x2,z | ^i] [z,x\ | Ifa] = [*2>*i I ^ = 
^1 U ^2] is ensured. This being the case, it is possible to associate each of the quark fields with its 
own gauge link because the attached contour has no bearing on the definition of fi/ a (x). Then, the 
struck quark can be replaced by an "eikonalized quark" 



W(x~\T) = \if(x~)[x~,°°~\r\ = v/t-O^exp 



-ig f X dz^(0 + ,Z-,0x)ta 

■J 00- [r] 



(1.3) 



which is a contour-dependent Mandelstam fermion field [Q] (with an analogous definition for the 
antifermion field). In this scheme, the gluon reconstitution in the gauge-invariant correlator for the 
integrated PDF involves gluons emanating either from the gauge links — giving rise to selfenergy- 
like diagrams — or contractions with the gluon self-fields of the Heisenberg operator for the struck 
quark which generate crosstalk-type diagrams. Note that for the sake of clarity and simplicity, we 
ignore bound states (spectators). As a result, one has 

£/fW = \Lj ^e- ik+ ^(P\%(^,0 ± \%) |«)7 + H^(0-,0J^ 2 ) |P) . (1.4) 
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This concept was carried over to unintegrated PDFs, i.e., to those PDFs which still depend on 
the transverse momenta — hence termed TMD PDFs. However, TMD PDFs with only longitudinal 
gauge links are not completely gauge invariant under different boundary conditions on the gluon 
propagator in the light-cone gauge A + = 0. The reason is that x~ -independent gauge transforma- 
tions are still possible under the same gauge condition. Hence, the naive collinear gauge-invariant 
TMD PDF definition as for the integrated case is inapplicable. Refurbishment is provided via the in- 
troduction of transverse gauge links which necessarily stretch out off the light cone to infinity [9, 0]. 
This generalizes Eq. ( |L~4| ) for a quark with k^ = (k + ,k~,k±) in a quark with = (p + , p~,0±) to 
the expression 

f qlq (xM = \j ^-^kexp(-^ + r+^-^) (q(p)W^-^ ± M-,^^-^±}' 



X [oo 



_] T 7 



,0 ± ;0-,0Jv(0-,0x)| 9 (p) 



(1.5) 



in which 



[oo-,£ ±;< r,£J = ^e X p 



[oo ,oo ± ;oo ,^ ± ] = ^exp 



P OO 

ig J dTn-A%t a (£,+rTT) 

ig f ' dT\-K a t a {£ L +\T) 
Jo 



(1.6) 



(1.7) 



are the lightlike and the transverse gauge link, respectively. 



2. One-Loop Gluon Virtual Corrections in the A + = Gauge 

The pursuit of a proper definition of TMD PDFs is a long-standing problem that was not ac- 
complished with the definition above. The reason is — frankly speaking — that nobody knows how 
the contour behaves at light-cone infinity when it ventures out in the transverse directions. This be- 
havior has influence on the singularity structure of the gluon propagator in the light-cone gauge 
A + = 0, notably, D^(q) = q i_^i +i0 (siiv — 9fi " v ^t| v " M ^ , via the boundary conditions to go around 
its singularities. To estimate this influence, one has to calculate the one-loop virtual corrections 
in the A + = gauge in conjunction with various boundary conditions (which absorb large-scale 
effects) and carry out the renormalization of the contour-dependent quark operators defined in Eq. 
(1.3). Two of us undertook this calculation, announced in ||,[8|], with a summary of the approach 



being given in [9]. The contributing diagrams are shown here in Fig. [T], while the corresponding al- 
gebraic expressions are given in Table [I] using the following symbolic abbreviations (the couplings 
g and g' below are labeled differently only in order to keep track of their origin; ultimately, they 
will be set equal): 

(i) Q: struck quark = e -Wdn¥^Wl y/f ee (£) — Heisenberg operator, 

(ii) longitudinal gauge link: [n~], 
(hi) transverse gauge link: [lj_], 

(iv) g refers to the QCD Lagrangian — see item (i), 

(v) coupling g' refers to the exponent of the gauge links, i.e., g' / °° dx ... , 

(vi) product g'g' corresponds to path-ordered line integrals in the exponent of the gauge links, i.e., 
g'g'IodzfQda... . 
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selfenergy crosstalk 

Figure 1: One-loop gluon virtual corrections to f q j q in the A + = gauge. The double lines describe the 
gauge links attached to the fermions (heavy lines), while the curly lines represent gluons, and the symbol <E> 
denotes a line integral. The Hermitian-conjugate (mirror) diagrams are not shown. 



Without going into too much detail, the results of this study show that the overlapping ultra- 
violet (UV) and rapidity divergences cannot be solely controlled by the dimensional (or any other) 
regularization. The ensuing divergence is of the type (1/e) ln(r/ / p + ), which becomes infinite when 
77 — > 0, and has, therefore, to be cured by an appropriate renormalization procedure. At this point it 
is important to mention that the terms on the diagonal in Table [I] represent selfenergy contributions, 
while all other terms are of the crosstalk type. In the gauge A + = only the terms QQ and Q[1J 
are non-vanishing. Moreover, the pole-prescription dependence in diagram (a) is canceled by its 
counterpart in (d) — see Fig. [TJ Taking into account the mirror contributions to (a) and (d) (not 
shown in Fig. [[]), one finds the following total contribution from virtual gluon corrections [||, [7p: 



Os,e) 



71 



1 3 , n 

- 7 +ln^- 
£ \4 p + 



y E + ln47T 



(2.1) 



struck quark longitudinal gauge link transverse gauge link 



struck quark 
longitudinal gauge link 
transverse gauge link 



QQ 



(a) 



Qk 



n n 



(b) 



Q[iJ 

[»-]PL] 



(d) 
(f)=0 

•(e) 



Table 1: Structure of the one-loop gluon virtual corrections to f q / q (x,k±) shown in Fig. 1. 
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From this expression one obtains for f q / q (x,kj_) the anomalous dimension ( y = \ 



T^-ioop = ^ C f f 7 + ln J = Smooth - Sy , (2.2) 



a s „ ( 3 , n 

where T] is the rapidity parameter with [tj] = [mass] and 8 y represents the deviation from the 
anomalous dimension of the gauge-invariant quark propagator in a covariant gauge (see [|3|] and 
earlier references cited therein). As argued in [2], [J], such an anomalous dimension can be 
associated with a cusp in the gauge contour at infinity and originates from the renormalization of 
the gluon interactions with this local contour obstruction. Therefore, one can claim that 8y can be 
identified with the universal cusp anomalous dimension |[I(]] at the one-loop order. But the choice 
of the gauge A + = should not affect the renormalization properties of the TMD PDF. Thus, the 



definition of f q / q (x,k±) given by Eq. ( |l.5| ) has to be modified by a soft factor (counter term) [ ]ll[ ] 

R = *(p+,n-\0)&(p + ,n-\£), (2.3) 

where <I> and are appropriate eikonal factors to be evaluated along a jackknifed contour off the 
light cone (the explicit expressions and a graphic illustration can be found in [|, g, fj). We have 
shown there by explicit calculation that in the A + = gauge with q~ -independent pole prescriptions 
(advanced, retarded, principal value), the anomalous dimension associated with this quantity ex- 
actly cancels 8 y, rendering the modified definition of the TMD PDF free from gauge artifacts. On 



the other hand, adopting instead a q -dependent pole prescription (Mandelstam [12], Leibbrandt 



130), no anomalous-dimension anomaly appears and the soft factor reduces benignly to unity 



3. Inclusion of Pauli Spin Interactions 



The conventional way to restore the gauge invariance of hadronic matrix elements is to use 



gauge links as those defined in Eqs. ( |1.6| ) and < \l.l\ ). However, this is only the minimal way to 
achieve this goal; it ignores the direct spin interactions because the gauge potential A^ is spin- 
blind. To accommodate the direct interaction of spinning particles with the gauge field, one has to 
take into account the so-called Pauli term ~ F^ lv S jlv , where S^ v = |[7^,7v] is the spin operator. 
Following this generalized conception of gauge invariance, we promote the definition of the TMD 



PDF to [14] 



fF/h(^ 



-Tr / dk' 



-*«<*|fc(«)[B-,«±;--,«_ 



xr[[. 



.;oo-,0 ± ]][[oo-,0 ± ;0-,0 ± ]] W -(0)|A) R , (3.1) 



where T denotes one or more y matrices in correspondence with the particular distribution in ques- 
tion, and the state \h) stands for the appropriate target. In the unpolarized case we have \h) = \h(P)), 
with P being the momentum of the initial hadron, whereas for a (transversely) polarized target the 
state is \h) = \h(P),S±). The enhanced lightlike and transverse gauge links (denoted by double 
square brackets) contain the Pauli term and are given, respectively, by the following expressions: 



[K,Ox;0-,Oj] = ^exp 



poo poo 

-ig do iiaA%(uo)t a -ig do Su V F^ v (uo)t a 
Jo Jo 



(3.2) 
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Symbols Expressions Figure 2 Value 

% J m dz\-£/{lz) (a) ^0 

^2 !qcIxS-^{ux) (b) ^0 

^3 / ~£/T5.Jf(lT) — 

^4 JodTjJdO (1-^(1(7)) — 

% JodzJ^da{\-^{\z)){S-^{\a)) — 

% J~dTj*do(S-3?(lz))(l-*f(lo)) - 

^ 7 ^dx£do (S-^(ux)) {S-^{uo)) (c) 

^ 8 £dT£da{S-&W) (5-^(1(7)) 

^9 f~dzJodo (W(1t)) {S-&{uo)) (d) ^0 

^io £ dz £ do (S-J^It)) (5-^(«a)) — 



Table 2: Individual virtual-gluon contributions appearing in the evaluation of Eq. ( 3.4 ) up to 



[[oo ,oo ± ;oo jOjJ] = ^exp 

3.1 Gauge links with Pauli terms up to 



/•oo /> oo 

ig dzl^AKizy-ig dzS^Frmf 

Jo ~~ Jo 



(3.3) 



Adopting this reasoning, we have to calculate in the A + = gauge the expression 

[[oo-,oo ± ;oo-,0j]-[[oo-,0 ± ;0-,0j] = l-/g(^ 1 +^ 2 + ^ 3 )-g 2 (^4 + ^5 + .-.^10) (3.4) 

with F a MV (oo-,fj+,£ ± ) = 0, Yi(^) = e^S^ ^tfyf*^) and an analogous expansion for the 
transverse gauge links (see Table ||), whereas the contributing diagrams are displayed in Fig. ||. Let 
us quote here some important features of the presented theoretical framework referring for details 



to our recent work in Ref. [14]: (i) The Pauli term is not reparameterization invariant — unlike the 



usual Dirac term. Therefore, we have to use the dimensionful vectors n+ — > m* = p~rit , — > 
Up = p + n~ , lj_ — > p + l±- (ii) The Pauli spin-interaction terms do not completely vanish along n~ 
in the A + = gauge, whereas terms containing JF(1t) (or J£"(la)) cancel out in the product of the 
gauge links and F fl ^ v (oo~,0 + ,£; ± ) = 0. (iii) To the g 2 -order level, the Pauli term reads 

5 • & = S MV J^ v = 2S+-^ + - + 2S +i ^ +i + 25_ l -^" _i + S t j& iJ (3.5) 

and has the following non-zero components: 

= d+^- ; = d+^i ; (3. 6 ) 

= d'srf 1 - d 1 ^- , & ij = d'a? j - d j £/> . (3.7) 

(iv) The diagrams (a)-(d) in Fig || represent virtual gluon corrections and contain UV and rapidity 
divergences that give rise to the anomalous dimension of the TMD PDF. In contrast, the diagrams 
(e)-(g), which describe real-gluon exchanges across the cut (vertical dashed line), contribute only 
finite terms. 

From Fig. ^ we see that the gauge-link correlator contains contributions of two different types 
related to selfenergy- and crosstalk-type diagrams. To discuss the structure of the correlator in a 
compact way, it is useful to use the following symbolic abbreviations: 
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Figure 2: One-loop gluon virtual corrections to f q / q in the A + = gauge. Graphs (a), (b), (c), and (d) 
describe virtual gluon corrections; graphs (e), (f), and (g) represent real-gluon exchanges across the cut 
(vertical dashed line). The double lines decorated with a ring represent enhanced gauge links containing the 
Pauli term. 

Q: Gauge self-field in the Heisenberg quark operator = v^Vl y/f ree (£) 

1 • stf (It) = A : Standard transverse gauge potential 

S ■ J^It) = F: Tensor (Pauli) term 

Then we obtain at (J{g~) the following results (consult Fig. || in conjunction with Table |2|): 
Selfenergy-type contributions 

• A X A X : (^ 4 )=0 not shown 

• F - F - :{%) = diagram (c) in Fig. g 

• F X F X : (%) = not shown 
Crosstalk-type contributions 

• QA X : (^} uv = -a,C F \ iC^ with C„ = {O(adv); -l(ret); -\ (PV)} — diagram (a). This 
term cancels the pole-prescription-dependent term in the UV-divergent part of the fermion 
selfenergy QQ. 

• QF~: (%) with (QF~)~ = (®f) and (QF - ) 1 - = {^) — diagram (b). Accordingly, for 
the leading twist-two TMD PDF, we find for the semi-inclusive DIS (SIDIS) 

r tw - 2 (%-) + (%-) f r^ 2 = ^c F r tw _ 2 , (3.8) 
rtw- 2 (^2 X ) + (^2 ± ) t r t w-2 = -^c F r tw _ 2 . (3.9) 

These two results combine to produce a constant phase (unrelated to that found in ^) 

St w - 2 = ct s C F 7t (3.10) 
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which is also valid for the twist-three TMD PDF, i.e., 6t w -3 = « s C F 7r, but flips sign for the 
Drell-Yan (DY) process because it depends on the direction of the longitudinal gauge link. 
Hence, our analysis [|l4|] predicts the important relation 

<Ssidis = -<5dy- (3.11) 

QF- 1 : (<%) = not shown 

A^F^: = —(^5), hence mutually canceling 

A^F^: (%)) = — diagram (d) in Fig. |2| — ("gluon mass" A 2 drops out at the end): 

( ^ 9) = "8^ CF[7+,r]r(£) ( 47r F) • (3 - 12) 
This nontrivial Dirac structure entails 



Tunpol. — 7 

Theiic. =7 + 7 5 

Ttrans. = Z^'V 



r unp oi. [y + ,y~] = -[7 + , 7~]r unpol . , (3.13) 
r h eiic. [7 + , 7"] = - [7 + , y"]r he uc. , (3. u) 

Ftrans. [Y + ,Y~] = -[Y + , Yl^ns , (3.15) 



where obvious acronyms have been used. Taking into account the mirror diagrams (not 
shown in Fig. Q), the twist-two terms mutually cancel by virtue of the relation 

[7+,7-]r tw _ 2 = -r tw _ 2 [7+,7-] = 2r tw _ 2 , 

which permits a probabilistic interpretation of the twist-two TMD PDF as a density on ac- 
count of A^F - — > 0. On the other hand, the twist-three TMD PDF gets a non-vanishing 
contribution to its anomalous dimension as one sees from 

r tw - 3 (%) + Wr^s = -J|[7 + ,7-]r(£) (a%^\ . 

• F-'-F - : (^10) = without assuming any particular form of the gauge field at light cone °°. 
3.2 Real-Gluon Contributions at G{ e 2 ) 



Besides the virtual gluon collections, there are also real gluon exchanges that contribute finite 
contributions to the TMD PDF. The main difference from the previously considered case is that now 
the discontinuity goes across the gluon propagator that has to be replaced by the cut one. More- 
over, the Dirac structures, marked above by the symbol T, are sandwiched between Dirac matrices 
stemming from Pauli terms standing on different sides of the cut. The real-gluon contributions are 
specified in Table 3. 

Using the same symbolic notation as in the previous subsection, we briefly remark that 

• F~F~: (^11) — > (at least power-suppressed ~ p~) 

. A^F-: (^i 2 ) + (^i 2 ) t ~r[7+,7-] + [7+,7-]r = 

• QF~: (^3) + (^3 y and (^3) + ) t mutually cancel up to a power-suppressed term. 
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Symbols Expressions Figure || 

%x Jo 'dz £ do {S-^{ux))Y{S-^{ua + ^-;^ A )) (e) 

% 2 fvdTf~do(l-^(lT))r(S-^(uo + t;-;S ± )) (f) 

jfc Sodor(S-&(uo + Z-;Z ± )) (g) 

Table 3: Individual real-gluon contributions to &(g 2 ) corresponding to the diagrams (e), (/), (g) in Fig. ^. 

4. Highlights and Conclusions 

We argued that the dimensional regularization of overlapping UV and rapidity divergences 
in TMD PDFs is not sufficient to render the TMD PDF finite — one needs renormalization [Q, 



To remedy this deficiency, a soft factor [jllj] along a jackknifed contour off the light cone 



was introduced into the definition of the TMD PDF [Q] whose anomalous dimension cancels in 
leading loop order the cusp anomalous dimension entailed by this overlapping divergence (with 
a full-fledged discussion being given in [^]). The modified TMD PDF reproduces the standard 
integrated PDF and is controlled by an evolution equation with the same anomalous dimension as 
one finds in covariant gauges with no dependence on the adopted pole prescription for the gluon 
propagator — this would be impossible without the soft renormalization factor (see [Q] and for a 
more dedicated discussion [|l"5|]). In particular, using the A + = gauge in conjunction with the 



Mandelstam-Leibbrandt pole prescription [12, AM, no anomalous-dimension defect appears and 



thus the soft factor becomes trivial. An important finding of this approach is that the anomalous 
dimension of the unpolarized TMD PDF for SIDIS and the DY process is the same, i.e., ySiDis = 
yPJ , albeit the sign of the e term in the gluon propagator ^+1^ is different for these two processes 
— irrespective of the boundary condition applied. Quite recently, Collins discussed alternative 



ways to redefine the TMD PDFs in such a way as to avoid rapidity divergences [|16[]. 

We also presented a new scheme for gauge-invariant TMD PDFs which includes the direct 
interaction of spinning particles with the gauge field by means of the Pauli term in the longitudinal 
and transverse gauge links. In some sense, the Pauli spin interaction is the abstract analogue of a 
Stern-Gerlach apparatus — sort of — and gives rise through the transverse gauge link to a constant 
phase 8 = a s C^%, which is the same for twist-two and twist-three TMD PDFs, but flips sign when 
the direction of the gauge link is reversed — thus breaking universality. As a result, one finds 
5dy = — ^sidis- To facilitate calculations, we developed in Ref. [ [14] ] Feynman rules for enhanced 
gauge links — longitudinal and transverse — which supplement those derived before for the stan- 
dard gauge links by Collins and Soper [|Xjl. Because the Pauli term contributes to the anomalous 
dimension of the twist-three TMD PDF, the evolution of such quantities is more delicate and may 
require the modification of the renormalization factor to preserve its density interpretation. 

Bottom line: Our results — most significant amongst them the appearance of a non-universal 
phase — may stimulate both theoretical and experimental activities. On the other hand, T-even 
and T-odd TMD PDFs may become "measurable" on the lattice, so that it seems possible that 
non-trivial Wilson lines, as those we discussed in this presentation, may be revealed in the future. 
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